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The magnetic and electronic properties of Ba4Ru30io were investigated by the ab-initio and model 
calculations. It is shown that nonmagnetic ground state of the one-third Ru 4+ ions is not due to the 
correlation effects. It is rather caused by the charge disproportionation between crystallographically 
different Ru and the molecular orbital formation in the Ru's trimer. 
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I. INTRODUCTION 

The complex Ru oxides are widely investigated in the 
last decades because of their unusual electronic and mag- 
netic properties. The odd- parity pairing was found in the 
superconductor Sr2Ru04, [l| while Ca2Ru04 became one 
of the model systems for study of the orbital-selective 
Mott transition. Q The integer spins on the Ru 4+ ions 
likely form the Haldane chains in TI2RU2O7, which leads 
to a drastic decrease of the magnetic susceptibility at low 
temperatures. [|| Similar behavior of susceptibility was 
recently observed in another Ru oxide - Ba4Ru30io- [3 
This drastic decrease was initially attributed to the for- 
mation of the unconventional antifcrromagnctic (AFM) 
state, but the mechanism of its stabilization is still un- 
known. 

Structurally Ba4Ru30io is made of the Ru-trimers. 
These trimers are formed by three RuOe octahedra shar- 
ing their faces. Each trimer is connected with four neigh- 
boring trimers (see Fig. [1} via corners of the outer RuC>6 
octahedra to build corrugated layers. We denote the mid- 
dle Ru ion in the trimer Ru m , while two outer Ru - Ru D . 

In Ba4Ru3 0io the Ru ions must have 4+ oxidation 
state with 4 electrons in the Ad— shell, since Ba is 2+ and 
O is 2-. Due to a large t^g — e g crystal-field splitting (see 
Sec. Ill for details) these four 4c?— electrons should be sta- 
bilized in the t^g sub-shell resulting in the state with S=l. 
Indeed, at room temperature Ba4Ru30io was found to be 
paramagnetic with the effective moment Heff = 2.83 [Lb 
typical for the S=l ions, [HQ while the decrease of the 
temperature below Tn = 105 K leads to the formation 
of the unconventional AFM state, where one-third of the 
Ru 4+ ions are not ordered (tentatively these are Ru m ). [H 

Several models were proposed in Ref. to explain this 
physical phenomenon. One may think that (1) due to 
some reason a part of the Ru ions stays paramagnetic 
down to a very low temperature; that (2) they are ran- 
domly frozen; or that (3) because of a strong spin-orbit 



coupling a nonmagnetic state of the Ru ion with J = 
is stabilized. 



In order to check all these possibilities, the ab-initio 
band structure calculations were performed for this com- 
pound. We found that the true reason for the absence of 
the magnetic moments on one-third of the Ru ions is the 
charge disproportionation between two crystallographi- 
cally noncquivalent Ru (Ru m and Ru D ) and formation 
of the molecular orbital in the trimers. In the band- 
structure picture it is manifested in the suppression of 
the nonmagnetic 4d density of states (DOS) of the mid- 
dle Ru m in the vicinity of the Fermi level. 
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FIG. 1: (color online). The crystal structure of Ba4RuaOio. 
Projection on the be— plane. Oxygen ions are in blue, Ruthe- 
nium - in red, Ba - not shown for the simplicity. Ru m are in 
the middle of the trimers, Ru - on the corners. 
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FIG. 2: (color online). The GGA partial density of states, 
obtained in the PWscf code. Ru m is the middle atom in the 
trimer, while Ru are the corner, outer to the trimers, ones. 
The Fermi energy corresponds to zero. 



II. CALCULATION DETAILS 

For the band structure calculations we primary used 
the pseudo-potential PWscf code. @ We utilized general- 
ized gradient approximation (GGA) with Pcrdcw-Burkc- 
Ernzerhof version of the exchange-correlation poten- 
tial @ and ultrasoft scalar-relativistic pseudo-potentials 
with nonlinear core correction (for better description of 
the magnetic interactions). The charge density and ki- 
netic energy cut-offs were taken 40 Ry and 200 Ry, re- 
spectively. 664 fc— points in a full Brillouin zone were 
used in the calculation. 

The generalized gradient approximation with the 
account of the on-site Coulomb repulsion (GGA+U 
method) was used to test an importance of the elec- 
tronic correlations. The Hubbard's on-site repulsion pa- 
rameter U and Hund's rule intra-atomic exchange Jh for 
Ru were calculated in Ref. H to be 3.0 eV and 0.7 eV. 

The Wannicr function projection was performed within 
the linearized muffin-tin orbitals (LMTO) method @ as 
it was described in Ref. [l^- 

All the calculations were performed for the crystal 
structure corresponding to T=10 K. Q 



III. NONMAGNETIC GGA CALCULATIONS 

We start with the analysis of the results obtained in the 
nonmagnetic generalized gradient approximation (GGA), 
presented in Fig. One may see that the oxygen 2p 
states are placed primarily from -7 to approximately - 
1 eV. The top of the valence band is formed by Ru-4d 
(t2 g ) states. The Ru-e g orbitals are above 1.5 eV. The 
^2g — £g crystal field splitting estimated by the center of 



FIG. 3: (color online). The band structure obtained in the 
nonmagnetic GGA calculations (in the PWscf code). The 
Fermi energy corresponds to zero. 



gravity calculation is 3.38 and 3.20 eV for Ru m and Ru Q 
respectively. Thus one would indeed expect that such a 
strong crystal field splitting would prevent occupation of 
the e g orbitals and d configuration of Ru 4+ would result 
in the S=l ground state of each Ru ion. 

However, there is another one quite important feature 
of the nonmagnetic DOS related to the splitting between 
the Ru m and Ru Q 4d states, which is clearly seen from 
Fig. [3J The same center of gravity calculation shows that 
this splitting exceeds 0.13 eV. The Ru m -4d states inte- 
grally turns out to be lower than Ru D -4d due to a larger 
mean Ru-0 distance. According to Ref. [3, the difference 
in the mean Ru-0 distances for these two Ru exceeds 
0.035 A. This is not a small number, since the difference 
in ionic radii between Ru 4+ and Ru 3+ is 0.06 A. [TO] As 
a result one could expect that there would be more d- 
electrons on the middle Ru, i.e. it would have on the 
average a smaller valence. And indeed, according to our 
nonmagnetic GGA calculations the difference in charge 
between Ru m and Ru D is 0.52 electron per atom, which 
is even larger than the charge disproportionation on the 
Fe sites in the charge ordered phase of Fe304. fl2l ] 

Thus, Ru m -i2 ff partial DOS states turns out to be al- 
most completely below the Fermi level, except a small 
shoulder and a narrow peak at ^0.9 eV. The detailed 
analysis of the band structure, presented in Fig. [31 con- 
firms that the band corresponding to this peak in DOS 
at 0.9 eV is dispersionless. One may also see that each 
band in the ZT direction is two times degenerate, since 
there are two formula units (f.u.) in the unit cell. In ef- 
fect each trimer provides one dispersionless band at ^0.9 
eV. 

In order to analyze the formation of this band we car- 
ried out the calculation within the local density approxi- 
mation (LDA) and the LMTO method for which one may 
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FIG. 4: (color online). The molecular orbital which cor- 
responds to the highest in energy Ru— ti g states (bands at 
~ 0.9 eV in Fig. [3}, obtained by the diagonalization of the 
three-site Hamiltonian with the use of the Wannier function 
procedure in the LDA LMTO calculations. 



perform the Wannier function projection procedure. [T(| 
The use of this procedure allows to get the on-site Hamil- 
tonian matrix for the Ru— Ad states and all the hoppings 
between different sites. The projection was performed 
onto thirty Ru— Ad orbitals (5 orbitals x 3 atoms x 2 
formula units) in the local coordinate system, where the 
axis are directed to the nearest oxygens. 

The diagonalization of the on-site Hamiltonian in a 
real space for each Ru leads to the order of the levels 
predicted by the local distortions. The tetragonal elon- 
gation of the Ru m C>6 octahedra results in two low-lying 
almost degenerate ^-levels split from the higher ti g or- 
bital by ^330 meV. However, a close inspection of the 
Hamiltonian shows that the hopping integrals between 
Ru G and Ru m exceed 360 meV, which is, in contrast to 
naive expectations, even larger than the off-diagonal on- 
site matrix elements for the t2 g orbitals. Thus, the or- 
bitals on which the Ru— Ad electrons tend to localize are 
determined not only by the crystal-field created by the 
local surroundings (oxygens), but also by the inter-site 
matrix elements. This is mainly due to a face-sharing 
packing of the RuOg octahedra in the trimcrs, which al- 
lows a direct overlap between d— orbitals on the neigh- 
boring Ru sites. Indeed, the nearest neighbors Ru-Ru 
distance across the common face 2.55 A is even less than 
that in ruthenium metal, 2.65 A. [l3| It is also quite im- 
portant that the relevant orbitals here are id— orbitals, 
which are more extended than 3d. fl4j 

In order to take into account the hoppings, we con- 
structed the large 15x15 Ru -Ru m -Ru Hamiltonian for 
the trimer in a real space. Diagonalizing this Hamilto- 
nian one obtains that the lowest and highest in energy are 
the molecular orbitals which have the a\ g symmetry with 
the largest contribution on the central Ru m ion. This is 
clearly seen in the Fig.[H where the highest in energy an- 
tibonding a* g molecular orbital is plotted. This orbital 
corresponds to the band, which has the energy ^0.9 eV 
in the ZT direction. 
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FIG. 5: (color online). The total and partial density of states 
(DOS) in the GGA calculation for the t m i — t m i mag- 
netic configuration (so called fully AFM), when the spins on 
the Ru both in the same and in the adjacent trimers are an- 
tiparallel. The results were obtained in the PWscf code. The 
Fermi energy corresponds to zero. 



Inspecting Fig. |U one may see that the contribution of 
the Ru m — Ad states to the antibonding a\ g orbital should 
be the largest (the same is actually correct also for the 
bonding a\ g wave function). This is easy to understand 
by solving the three-site (i = 1,2,3) problem with the 
single orbital per site in the tight-binding approximation. 
For simplicity the energies of each orbital are taken the 
same, Si — 0, and the hoppings, t > 0, are nonzero only 
between the nearest sites. For the linear cluster of three 
atoms the Hamiltonian is written as 




(1) 



The energy spectrum consists of the bonding (b), anti- 
bonding (ab) and nonbonding (nb) orbitals: 

E ab = tV2, ibab = -(ipi - V2> 2 + V's), (2) 



E nb = 0, 1p nb = — ^(V>1 ~ -03), 



(3) 



E b = -ty/2, 4> b = + y/2fo + (4) 



As a result one may see that the weight of the states 
corresponding to the middle, i = 2, atom to the bond- 
ing and antibonding wave functions equals to 50%, while 
they do not contribute to the nonbonding molecular or- 
bital at all. Thus, the Ru m — ti g states of the a\ g sym- 
metry are pushed away from the Fermi level due to the 
molecular orbital formation. Moreover, the shift of the 



rest of the Ru, 



tig states downwards, with the cor- 
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TABLE I: Total energies, magnetic moments on the Ru and 
Ru m ions and band gaps. Results of the GGA in the PWscf 
code. Total energies are per formula unit (f.u.). The notation 
f m 4 means that two Ru ions are AFM, while moment for 
the Ru m can be found in the 3rd column. Symbol "-" denotes 
the bond with adjacent trimers. The unit cell vectors are cho- 
sen in such a way, that the type of inter-trimer magnetic or- 
der (ferro- or antiferromagnetic) for presented configurations 
is the same for all 4 inter-trimer bonds. 





Total energy 
(meV) 


Spin moments 
Ru D /Ru m (/i B ) 


Band 
gap (eV) 


t m 4 


- 1 m 4 





0.95/0 


0.12 


t m t 


-tmt 


22.7 


0.98/0.33 


Metal 


t m 4 


- 4 m t 


20.1 


0.94/0 


Metal 


fmt 


— 4 m 4 


61.0 


0.69/0.15 


Metal 


OmO 


- OmO 


137.6 


0/0 


Metal 



responding charge disproportionation, leads to the sup- 
pression of the Ru m Ad partial DOS in the vicinity of the 
Fermi level. The Ru D — 4c? states in contrast form well 
defined peak exactly at the Fermi energy. 

The Stoner criteria cannot be applied directly to the 
partial DOS, but it is clear that the gain in the magnetic 
energy due to the spin splitting on Ru m will be minimal 
(because of the pseudogap in its partial DOS), while the 
loss of the band energy will be substantial. Thus in the 
band picture the Ru m — Ad states are not expected to 
be polarized. In contrast, the Ru — Ad states can be 
easily magnetized. This is exactly what we observe in 
the magnetic GGA calculations. 

IV. MAGNETIC GGA CALCULATIONS 

Several configurations were investigated in our mag- 
netic GGA calculations (see Tab. HJ. The lowest in en- 
ergy is the state when the spins on the Ru D both in 
the same and in the adjacent trimers are antiparallel 
(f to I — t m I)- The partial DOS corresponding to 
this configuration are plotted in Fig. [5] One may see 
that there is no splitting between majority and minority 
spins of the Ru m . Thus according to our calculation 1/3 
of the Ru atoms indeed are nonmagnetic, exactly as it 
was observed in the experiment. [3] 

The magnetic moment on the Ru m appears only in 
the configurations with the FM order of the Ru D in the 
trimers. The ferromagnetically ordered spins of the Ru G 
create an exchange field on the Ru m , which slightly mag- 
netize these ions. On the band-structure language the 
appearance of some magnetization on Ru m in a ferro- 
magnetic state can be explained by the spin splitting of 
the Ru m — Ad states due to their hybridization with the 
Ru G — Ad states, which as we have shown above is large 
enough. 

The magnetic moments on Ru Q are smaller than 2/j,b 
expected for S=l. This is a rather typical situation for 



the transition metal oxides (especially for those based 
on the 4d and 5d elements) and can be attributed to a 
strong hybridization with the oxygens. Because of the 
hybridization we also observe non zero magnetic mo- 
ments on some of the oxygens, which vary from 0.04 to 
0.26 hb, depending on the magnetic configuration and 
oxygens crystallographic positions. The largest magnetic 
moment, 0.26 fis, was found on the oxygen atoms, which 
belongs to two neighboring trimers in the fully FM solu- 
tion (t m t - t Tn t). 

The analysis of the calculation results presented in 
Tab. U shows that the magnetic contributions to the total 
energies of different solutions is not described by purely 
Heisenberg terms, like JijSiSj. Indeed, if one compares 
the total energies of two solutions with the AFM order in 
the trimer (e.g. f m 4- — t m i an d t 171 1 — I m t) then 
the lowest will be the one with the AFM coupling be- 
tween trimers. However, if one analyses another two solu- 
tions with the FM order in the trimer (e.g. f m t — t m t 
and t to t - 4 m 4) , it becomes clear that the lowest will 
be the one with the FM coupling between trimers. This 
is mostly related with the smaller magnetic moment on 
the Ru m ions in the j" to j" — 4 rn 4 configuration, and it 
may be a consequence of molecular orbital formation on 
trimers" . 

There are also two other contributions to the magnetic 
energy in addition to the conventional Heisenberg model, 
which have to be taken into account to describe magnetic 
properties of Ba4Ru30iO. This is the development of the 
magnetic moments on Ru m (when the spins in the trimer 
are parallel) and on those oxygens, which are shared by 
two trimers (when the spins on the neighbor trimers are 
FM ordered). In each case the contribution to the total 
energy is —IM 2 /A, where M is the magnetic moment on 
the Ru rra or O ions, and / - the Stoner parameter, which 
is approximately equals to the Hund's rule exchange pa- 
rameter for Ru (0.7 eV) and to 1.6 eV, as it was calculated 
in Ref. pj] 

It is important to mention that only fully AFM state, 
corresponding to the lowest total energy, is insulating in 
the GGA approach, without including electronic corre- 
lations (Hubbard's U). This is related to the band nar- 
rowing caused by the antiferromagnetism. The small gap 
insulating ground state, obtained in the magnetic GGA 
calculations agrees with a semiconducting temperature 
dependence of the electric resistivity, observed in Ref. Q 



V. CORRELATION EFFECTS 

It is well known that for the description of the elec- 
tronic and magnetic properties of the transition metal 
oxides one often needs to take into account strong 
Coulomb correlations, which can be incorporated in 
the calculation scheme via the LDA+U/GGA+U Q or 
LDA+DMFT [3 formalism. 

We repeated the total energy calculation for some of 
the magnetic configurations within the GGA+U approx- 
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TABLE II: The total energies of the different magnetic solu- 
tions, as a result of the GGA+U calculations in the PWscf 
code. Total energies are given per formula unit. The nota- 
tions are the same as in Tab. [T] 





Total energy 
(meV) 


Band gap 
(eV) 




- 1 m 4- 





0.26 eV 


t m t 


-tmt 


18.3 


Metal 


t m I 


- jmt 


19.8 


0.09 eV 


t m t 


- |m j 


86.4 


Metal 


OmO 


- OmO 


284.6 


Metal 



imation (as mentioned above, we took the values U=3 eV 
and Jjj=0.7 eV, obtained in Ref. 0) and found that the 
account of the Coulomb correlations increases the band 
gap values, but does not basically change the order of the 
total energies for different magnetic configurations (com- 
pare Tab. HI and HI)) . This confirms that the origin of the 
unconventional magnetic properties is not related to the 
correlation effects, typical for the localized electrons, but 
is largely due to two factors discussed above: the for- 
mation of molecular orbitals on Ru trimers and charge 
redistribution, with the increase of d-electron number on 
the middle ruthenium, Ru m . 

VI. SUMMARY 

Summarizing, on the basis of ab-initio and model cal- 
culations we obtained the explanation of the unusual 



magnetic properties of Ba4Ru30io, which consist of the 
Ru trimers coupled via corner-shared oxygens. The 
ground state is found to be an antiferromagnetic insu- 
lator, with antiferromagnetic ordering both within and 
between Ru trimers. The most surprising fact - the non- 
magnetic nature of the middle Ru ions in each trimer - 
is explained to be a result of a combined action of the 
formation of molecular orbitals in the Ru trimers and of 
charge redistribution, with the increasing d— states occu- 
pation on the middle Ru which is sandwiched between 
two other Ru ions with antiparallel spins. These unusual 
magnetic properties of Ba4Ru3 0io are well explained by 
the band picture, and the electron correlations do not 
play significant role in the formation of this magnetic 
state. 



VII. ACKNOWLEDGMENTS 



We thank Igor Mazin and Zlata Pchclkina for the 
fruitful discussions about band magnetism. This work 
is supported by the Russian Foundation for Basic Re- 
search via RFFI-10-02-96011 and RFFI-10-02-00140, by 
the Ural branch of Russian Academy of Science through 
the young-scientist program, the Ministry of education 
and science of Russia (grant 12.740.11.0026), by the Ger- 
man projects SFB 608, DFG GR 1484/2-1, FOR 1346, 
and by the European network SOPRANO. 



[1] K. D. Nelson, Z. Q. Mao, Y. Maeno, and Y. Liu, Science 
(Ne w York, N.Y.) 306, 1151 (2004), ISSN 1095-9203 
URL |http : //www . ncbi . nlm . nih . gov/pubmed/ 1 5539595 [8] 

[2] E. Gorelov, M. Karolak, T. O. Wehling, F. Lecher- 
mann, A. I. Lichtenstein, and E. Pavarini, Phys. 
Rev. Lett. 104, 225401 (2010), ISSN 0031-9007, URL [9 ] 
http : //link . aps . org/doi/10 . 1103/PhysRevLett . 104 . 226401 

[3] S. Lee, J.-G. Park" TT Adroja, ET Khomskii, 

S. Streltsov, K. A. McEwen, H. Sakai, K. Yoshimura, [10] 
V. I. Anisimov, D. Mori, et al., Nature ma- 
terials 5, 471 (2006), ISSN 1476-1122, URL 
|http : //www . ncbi . nlm . nih . gov/pubmed/ 16699512 

[4] Y. Klein, G. Rousse, F. Damay, F. Porcher, [11] 
G. Andre, and I. Terasaki, Physical Review 
B 84, 054439 (2011), ISSN 1098-0121, URL 



Rev. Lett. 77, 3865 (1996), ISSN 1079-7114, URL 
http://www.ncbi .nlm. nih. gov/pubmed/10062328 
V. Anisimov, F. Aryasetiawan, and A. Lichtenstein, 
J. Phys.: Condens. Matter 9, 767 (1997), URL 
http : //iopscience . iop . or g/0953- 8984/9/4/0(?2| 
O. K. Andersen and O. Jepsen, Phys. 
Rev. Lett. 53, 2571 (1984), URL 
|http : //link . aps . org/doi/10 . 110 3/PhysRe vLett . 53 ■ 2571 | 
S. V. Streltsov, A. S. Mylnikova, A. O. Shorikov, 
Z. V. Pchelkina, D. I. Khomskii, and V. I. Anisimov, 
Phys. Rev. B 71, 245114 (2005), ISSN 1098-0121, URL 
|http://linkTa ps . org/d oi/10 . 11 03/PhysR evB . 71 ■ 245114| 
R. D. Shannon, Acta Crystallograph- 

ica Section A 32, 751 (1976), URL 
|http : //dx . doi . org/10 . 1107/S0567739476001551 | 



http : //link . aps . org/doi/10 . 1103/PhysRevB . 84 . 054439 

[5] C. Dussarrat, F. Grasset, R. Bontchevb, and J. Darriet, 
Journal of Alloys and Compounds 233, 15 (1996). 

[6] P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, 
C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococ- 
cioni, I. Dabo, et al., J. Phys.: Condens. Matter 21, 
395502 (2009). 

[7] J. Perdew, K. Burke, and M. Ernzerhof, Phys. 



[12] L Leonov, X. 1NL Yaresko, V R Antonov, 
M. A. Korotin, and V. I. Anisimov, 
Phys. Rev. Lett. 93, 146404 (2004), URL 
|http : //link . aps ■ org/doi/10 . 1 103/PhysRevL ett . 93 . 1 46404 

[13] E. A. Owen and E. W. Roberts, Zeitschrift fuer Kristal- 
lographie, Kristallgeometrie, Kristallphysik, Kristall- 
chemie 96, 497 (1937). 

[14] J. B. Goodenough, Magnetism and the Chemical Bond 



6 



(Interscience- Wiley, 1963). 
[15] I. I. Mazin and D. J. Singh, Physi- 
cal Review B 56, 2556 (1997), UR L 
|http : //prb . aps . org/abstract/PRB/v56/i5/p2556_l | 



[16] V. I. Anisimov, A. I. Poteryaev, M. A. Korotin, A. O. 
Anokhin, and G. Kotliar, J. Phys. Cond. Matt. 9, 7359 
(1997). 



